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ABSTRACT 


The model of a Surface Effect Ship was refined to include 
Simplified propulsion dynamics, negative drag characteris- 
MmiGoye sca Stdee CLLrects and an autopilot for speed control. 
These design modifications were introduced into a real time, 
ome ontrolled Simulation of a 3000 ton Surface Effect Ship 
(3K-SES) in 5 degrees of freedom (RTS5D) and results were 
compared with a Data Base Program (DBSIMS5D) based on towing 
tank data scaled up to model the ship. 

Hardware and Software design changes were incorporated 
into the RTSSD model to provide a more accurate approxima- 
Euameot real time, a faster computer iteration time, and a 
meagre condition Warning 1f the operator exceeded certain 


mimGist Vectoring limits. 
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ae) CRODUGHION 


Besime tuel Costs, Critical manpower shortages and an 
increasing need for rapid deployment of forces in defense 
of United States' interests throughout the world mandate 
the development of fuel efficient, high speed, minimum man- 
ned ships for the Navy. The Captured Air Buddle (CAB) Sur- 
Eecmercet ship (SES) holds great promise in providing an 
answer to these problems. 

A ship with minimum manning traveling at the high speeds 
emvrsioned for the SES would require highly skilled operator 
personnel. It would be necessary to provide these personnel 
with real time operational training similar to that which 
merots receive in aircraft simulators. This would test the 
man in responding to failure situations that may be encoun- 
tered during actual operation and train personnel in control 
techniques that are unique to SES type vehicles. A Real Time 
Simulator could also be used as a means to test hardware used 
mimoko Cratt and provide a design tool for the development 
of future modifications. 

A simplified, non-linear five degree of freedom Real Time 
Simulator (RTSSD) for the 3K-SES has been developed by T. S. 
Newson (Ret. 1] such that a continuously observable real time 
solution is interfaced with man-generated control. Certain 


'man-in-the-loop'" experiments have been conducted using the 


Acs 





RTSSD that have demonstrated limitations in that model. It 
was necessary to: 1) refine ship dynamics; 2) introduce sim- 
plified propulsion dynamics; 3) constrain thrust angle in 
order to remain within the bounds of programmed dynamic res- 
ponse (avoid broaching); 4) refine the software to provide 
faster computer iteration time; and 5) change the hardware 
to provide the operator with speed control as well as turning 
Sontro]. 

mecomaingly, changes have been made in the RTSSD model 
mmiore dectinately reflect the nonlinear forward drag 
characteristics of the SES. A speed controller has been 
designed to allow operation of the 3K-SES within a speed 
range of 40-60 knots and limits have been placed on the 
thruster deflection angle in order to prevent broaching, a 
condition which is not allowed under the existing dynamic 
model. Software refinements have been incorporated to im- 


prove computer iteration time. 


» 
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Il. EQUATIONS OF MOTION 


A. COORDINATE SYSTEMS AND ASSUMPTIONS 
The RTSSD model of the 3K-SES used simplified equations of 
motion. For completeness of comparison the development of 
these equations is repeated here. The rational for requiring 
Simplification was to start with differential equations of 
motion that could be utilized effectively in a real time simu- 
lation. To satisfy this requirement a model was selected that 
would use point source of drag forces consisting of only a 
few lumped-coefficient terms. The equations used are based 
upon those developed in Ref. 2 which follow this concept. 
Additional assumptions used in the simplified 5 DOF equations 
were as follows: 
Mmextimedeceleratr1ons arewWeasured at the center of 
gravity. 
2) Vehicle is "free-to-heave". 
Smetlimeress COUDLed MOMents Of iInewela are Zero: 
ee ewermoss progmets of angular velocities in force 
Squatrons are Zero. 
5) All roll and pitch angles used in force equations 
are subject to small angle approximations. 
6) Calm water conditions. 
Mie ccc Ndassecistrtpution OL vehicle are Constant. 


8) Effect of changes of aerodynamic force are neglected. 


LS 





Utilizing the above assumptions, the simplified 5 DOF 


equations of motion are developed using the coordinate sys- 


tems shown in Fig. 1 and Fig. 2, the force vector diagrams 


Shown in Fig. 3, Fig. 4A, Fig. 4B, and the following 


equations from Newton's laws of motion: 


SURGE 
SWAY 
YAW 
EoD er 


ROLL 


where: 


m 


m(u - vr) 


ae - ur) 


i Ml 
< ON 


I 
Pa 


= M 


= K 


mass of the rigid ship 


Velocity elit x - aI recetom (SURGE) 


velocity in y-direction (SWAY) 


aieelar Velocity about they z-ac1s 


angular velocity about the x-axis 


angular velocity about the y-axis 


moment of 
moment of 
moment of 
summation 
summation 


summation 
Z=axX1S 


summation 


summation 


inertia about the z-axis 
inertia about the y-axis 
inertia about the x-axis 
of forces in x-direction 
of forces in y-direction 


of moments about the 


of moments about the y-axis 


of moments about the x-axis 
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eee seas 
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ft-1b. 


ft-lb, 


ft-lb, 





Additionally the following navigation relationships are 


Miaiazed {see Fig. 1). 


1 =u cosp - v Siny 
Y5 aa Sing + Vv cost 
V oe a : a 
S 
a ae 
8 = tan (-v/u) 
where: 


vy = heading angle 


ees Curt angle 


6 = thrust vector angle 
ESO , 

X = rx dt = ey ede 
O O 

x = sY dt mera 
O O 

iP = aoa cd Ve) mie Ce 

p = fp dt Ose Om b 


ee 








O TNAV 


Note: 1) w=0 when vehicle heads paralled to + Xray 
axe S. . 


2) Vehicle in figure above shown with negative 
thrust vector angle 6, positive wv, negative 
pppe Zative sway Velocity Vo01-e-, an a nicht 
turn 


Pelee: 


Definition Of Coordinate System (Part I) 
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Note: Direction for positive yelocity and 
angular rates are shown 


alk Wem 2 


Definition Of Coordinate System (Part ITI) 


LE, 





HCyyv iv | 





STARBOARD 
3 IDEWALL 


PORT 
SIDEWALL 


£ T5217 15 


Note: 1) p and p' are equivalent point force centroids. 
2) So; 5.» Sg» S56 Bre wChEbsite weetor anoles ang 
are not required to be equal. 
5) TA, Tg 5 Tg» Ty are thrust magnitudes and are 
not required to be equal. 
4) Force directions shown are for a right turn 


Figure 3 


Surface Effect Ship (Top View) 


20 








Sece: 1) Right Turn Forces Acting, i.e., Sa 5g > Sg: 
$49 are negative 


2) Torne = UPDRS coal Ue ee Oe Nae a 


Figure 4A 


Surface Effect Ship (Stern View) 
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Note: T =Poeso. + l--cOS6.,+ i xeoso.+ !, .coss 


FORW ~7 i aes § 9 9 “10 
(see Fig. 3) 
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Rouge .4.5 


Surface Effect Ship (Side View) 
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B. FORCES AND MOMENTS 
1. Surge Forces (see Fig. 4B) 

The forward acceleration was determined to be genera- 
ted by the summation of thrust vectored in the forward direc- 
tion ore The resultant acceleration of the SES was 
assumed to be opposed by a retarding force exhibiting a 
‘velocity-squared" characteristic. Additionally, Newton's 
laws of motion specified a '"'v.r'' product that contributed 
to the drag component of the equation when the vehicle was in 
a turn (contrifugal force reaction term). The simplified 
equation of motion for the surge acceleration follows from 
a summation of forces in the u direction as depicted in 


Ee, 48, 


o 


_ (Te. /m) - (Cpyu"/m) ge 


2. Sway Forces (see Fig. 4A) 
The side acceleration of the vehicle was analyzed 
Eee pe the result of the summation of the thrust vectored 
parallel to the vehicle's y-axis and a retardation "'velocity 
squared'' phenomena such as that produced by the cross flow 
drag term of the sidewall as described in Ref. 3. This 
retardation force is augmented by a "u-r" product as speci- 


fied by Newton's Law. 


Vv = Tosg_/m - Cpy vivi/m - ur 





Note the requirement to force the Cayo /m term to 
maintain the sign of the sway velocity such that both left 
and right turns may be computed. 

53. Yaw Moments (see Fig. 3) 

The yaw acceleration was a turning moment summation 
which was defined to be influenced by three primary forces 
Pperating Over moment arms. It is of interest to note that 
the third term is an "added mass" term which is considered 


to operate on the same pressure point as the C,.v|v| term. 


Di 
The term Taw 1S a Summation of moments generated by the 
PeollGethrusters as shown in Fig. 3. 
aw = S4T4 cosé. + S31, COS 5g - $417, COS Sq 


S5T 15 cos 540 - aot. Sin 55 - 2 ots sin 5 


- 21, sin Sg - aot i0 sin 819 


= Chy ft vivi/t, + A,juv2 /I, + T 


DY 22 yaw! 12 


4. Pitch Moments (see Fig. 4B) 
The pitch acceleration was assumed to be the summa- 


tion of moments generated by forward thrust, the 


Veorw? 
buoyancy of the sidewalls of the SES, Fy and Poe and a 
vertical force generated at the bow of the vehicle, F.. 
This vertical bow force was defined to be a lumped parameter 


term which modeled the reaction force due to plenum pressure 


acting against the bow seal. 
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Maan YD i 3 


average draft of bow seal 


z 
x2 + (C2, tand)/2)2, 


i 
"xx, 2 (2, tane)/22_ 
A e pg 
i 2 
A Q pg 
> > 
Eee h ed 31 ot 


average draft of sidewall 
plenum pressure 


average weteed sidewall areagof the stern 
width of one sidewall 


width of bow seal 


average wetted sidewall area of the bow 


agdea Mass COCErICiciEe 
water density 

S255 ft/s* 

lever arm of bow seal 


lever arm of buoyancy force 


p23 ~ Azgual/Ty 


ZS 





Note the added mass term A, zuq which was required 
to provide damping to the pitch moment and is specified 
mm Ret, 2. 

It was found that the response of the SES to a 
pitch perturbation without the added mass component was 
undamped and approximately sinusoidal. Unlike the yaw 
equation where the added mass term was a small contributor 
to damping, the pitch added mass term was found to be of 
Bieniticant importance in modeling the known pitch motion. 

2. Roll Moments (see Fig. 4A) 

The roll acceleration equation was analyzed to be a 
Summation of moments generated by buoyancy forces of the 
port and starboard sidewalls, ao and Bes the thrust 
wmeetared parallel to the y-axis Ometne= Vehicle, Tes a side 
force Cp. apace wawemiped oer ictomeuvercical 1oree, 

P we mirenvcreical force Eee was defined as the force 


generated in a turn due to the sidewall curvature (dead 


[meen anglé) acting against the cross flow of water. 


= = pe Se ae 
Ss ate ae 
7 cpze * |! 
2a aaveraoe drdt tvotuskouestacwall 
it 
a a, ee 
= oh tj 
* ds am on oY 
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te = average wetted area port 


ws average wetted area starboard 


emt ss) sw side ay “DY viv] 2, 


: Sew ae ~ Agy up)/Ty 


Again note the added mass term A.zup/Ty which was 
found to be essential in modeling the damping phenomena of 
eaes roll motion. 

In summary, the simplified 5 DOF equations of motion 


Meedg in RISSD for the 3K TON SES are: 





SURGE u = eee - Cy ae 1 oer Qe rey 28) ) 
SWAY v= Tide/™ - Cny viv{/m = ur (Fig. 4A) 

YAW t= Tyay/lz + Cpy vivi 2/1, + Ayguve,/T, 
Geis 3S 

PITCH a = (Teorwtp * F5ts1 + Fots - Spy lee 

- Fi 2. : Az,uq)/Ty (Fig. 4B) 

noe p " (CF. eee ose TSide*p * Cpy vivit, 
SF up)/Iy (Fig. 4A) 
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Peer ARAMETER IDENTIFICATION 
These equations are an extension of the 3 DOF flat turn 
SES model developed by Gerba and Thaler in Ref. 4. The 


identificati , 
tion of craft parameters Chy Chy> and 2 is 


described in Ref. 4 and repeated here for completeness. 


The surge drag coefficient C is determined by selecting 


DX 


a steady state turn condition, where T Mee eG eae ont 


Pog Ws 
are known. 


Z 
a Goel 
, forw DK, & 
m m 
mom which 
: _ Ss i aa ales 
DX 


u 


ite sway drag coefficient Chy is determined by the same 


steady state turn condition which requires 


: Tes C 
Weer = —2Se 2 DY yiyi - ur 
™m m 
where 
C side pee 
DY Viv 


The sway drag moment arm follows utilizing the yaw accel- 


eration equation 


Q 
Pe ae . Coyvivl2,, ; Ags uve, 
I, 1, I, 


28 





which yields 


-T 
yaw 


eeeeviy| + A 


DY DI 


The addition of the pitch and roll equations introduces 


additional lumped coefficients Cyp and Chzp: These are 


easily solved using known constants 2 Q Q 
De 5 Ww? oe Chy> z? dee 


We» Phe 2g and steady state values of eee Hoes Fy» F,; 
ea ew ? ee ee eee, and Tf. Lt IS tsionitteantes co eno Fe wena 
ima steady state condition the angular roll rate, p, and 
Beeular pitch rate, q, are both required to equal zero. 
This is in contrast to the yaw rate, r, where a finite 
Steady state value is desired-in a turn. Thus the added 
[yD uN eles jolie la ete me) ly eee) — 


mass terms of A 3uq and A 


3 om 
tions are not utilized in the determination of craft para- 
Meeers;, their function 1s strictly confined to damping of 


@itelr respective accelerations. Therefore, utilizing the 


pitch equation with known steady state turn values. 
od = (7 h + CopPpWe (2g - %3, tan 9)25, 
Ea ey a 
2°35 DX lee 3 af 


yields 
C ws be 2 Siar ee 


C eX eS forw * 
DP PW (24 = hey tan B)224 


c 


z9 





The roll equation under steady state turn conditions 
yields 

— 0) — - { = bs 

Dp UL Fe, Bee) ee eae Cay viv| 22 


- Coop viv et) 1 = 


from which 


“Dz ; CLES» : ee a TSide*p CDy Mary ge eo 
ee ee ene 
where 

F op = port sidewall buoyancy in pounds (Fig. 4A) 

Ec = starboard sidewall buoyancy in pounds (Fig. 4A) 
ee = thrust vectored parallel to vehicles yee 

in pounds fees 7 oy 

ie = sway steady state velocity in feet/sec 

Cny seoWay ocdrac coetitcient for steady state Conditren 
e = sway drag moment in feet (Bac 
- = effector thrust moment arm in feet (Fig. 4A) 
a, = surge drag moment arm in feet (Fig. 4A) 
Q = sidewall buoyancy force moment arm 

7 in feet (Fig. 4A) 
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A. DBSIMSD BENCHMARK 

Validation of the RTSS5D was attempted using the data 
base program (DBSIMSD) [Ref. 5S] as a benchmark. However, 
new guidance in the proper use of the data base program con- 
tained in Ref. 6 revealed the following errors in the use of 
the DBSIM5D. 1) Initial conditions had been improperly set 
and 2) Data used for validation for the 56 knot run at 
eeec@eor deflection angles greater than 10° was not correct 
because of the thrust inlet broach conditions which existed 
using the DBSIMSD. These conditions were not accounted for 
imcie RISSD model. The noted discrepancies in the opera- 
tion of the DBSIMSD necessitated a revalidation of the 


RiSoD. 


PeeevaLiLNATITON RESULTS 

isenme ene information contained in Ref. 6, the DBoIM5D 
program was again used as a benchmark to validate the RTS5D 
simulation. A DBSIMSD sequence of runs with initial forward 
velocities equal to 40 knots, 50 knots and 60 knots with 
meoomettector angles of 5° and 10° (at 40 knots), 5°, 10° 
amas: (at SO knots), and 5° (at 60 knots) were compared 
mumeideneical runs of the RISSD. Thrust effector angles 


were kept below a maximum that would have caused a broach 


Dal 





condition. Broach conditions are explained in Chapter IV. 
The first peak overshoot value and "Quasi Steady State" 
values for the variables u (surge), v (sway), r (yaw rate), 
¢ (roll), and 6 (pitch) were used as a basis for comparison. 
mivesresults are shown in Tables I-III. Significant error 
in forward velocity occurred in the 50 knot test. Also, it 
was noted that in 40 knot test with 5° thruster deflection 
angle that the DBSIM5D actually reached a ''Quasi Steady 
State" forward velocity greater than the initial value upon 
entering the turn. These conditions led to the development 
of the model of the RTSSD with negative drag characteristics 
in sea state. Differences in the other measured parameters 
were addressed by parameter adjustment on the new model in 


order to more closely match values obtained with the DBSIMSD. 


or 





TAREE ol 


RTSSD and DBSIM5D Performance Test at 40 Knots 


effector 
angle DBSIMSD RTSSD 
> in degrees o 
lst pk Qss lst pk Qss 
a er) ee Le le 
£0 ° 1.46 roo ol: 24 
> in degrees > 
lst pk Qss lst pk Qss 
a 1.44 ai a lee 
PO ° 1.48 1.48 eaves Bea 
u in u 
lst pk Qss is te poms Qss 
5 67.4 70.9 n/a 67.03 
1G ° my a Bio ny a 65.98 
V in V 
Est pk Qss ist apk Qss 
a 1.42 eae 2.84 Cae, 
> oe, 2.0 Se are o 3.41 
a in degrees/sec x 
ese pK Qss eres pik Qss 
a . 396 ~403 . 58 a 
10 955 959 .98 66 
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TABLE [I 


RTS5D and DBSIMS5D Perfermance Test at 50 Knots 


effector 
angle DBSIMS5D Rae) 
> Wy le wires > 
lst pk Qss lst pk Qss 
5 ~424 .446 ie eg 
0 ° La . 886 49 J o8 
HES ~ ea ie a oe 
6 in degrees 6 
St pk Qss lst pk Qss 
5 S02 .892 ila ie HS 
O° .943 nies eal,.5 Ihe Le: 
5 ieee eS Lee bey 
u Ite mee see u 
St pk Qss Es Gypk Qss 
5 ay a S540 my a 83.8 
eg) ° n/a oe n/a S2aro 
feS ° n/a 693.2 n/a 80.58 
V ii ekey sec V 
lst pk Qss PS pk Qss 
5 ie) ioe 5559 55 GZ 
eo ~ el 2 2.70 4.79 4.26 
5° 4.05 a So bo eS See Oy 
AG in degrees/sec i 
lst pk Oss lst pk Qss 
Die . 368 DEI! aS a 
EO ~ . 864 ow) eee oe 
#5° Lis 1.44 OG P26 
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TABLE ii! 


RTSSD and DBSIM5D Performance Test at 60 Knots 





eerector 
angle DBSIMSD Riso 
p in degrees > 
lst pk Qss PS Camis Qss 
B 141 2200 Bog p74 
8 in degrees 8 
lst pk Qss Ist pk Qss 
a oe 00 1.14 ie 7 
u in ft/sec u 
Fst pk Qss lst pk Qss 
os n/a 97.8 n/a 100.29 
V in ertysec V 
at. pk Qss lst pk Qss 
oe 205 1.84 4.26 oO 
r in degrees/sec iG 
lst pk Qss lst pk Qss 
oa Poo8 .494 Oe .487 
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IV. MODEL DEVELOPMENT 


fee DRAG CHARACTERISTICS 
1. Surge Equation Model 

The RTSSD model assumed a velocity squared drag re- 
tarding characteristic for the 3K-SES which led to the equa- 
tion for surge acceleration that was presented in Chapter II. 
iesdrag characteristics from Ref. 7 for the 3K-SES that is 
Simulated by the DBSIMSD are shown in Fig. 5. The data base 
program was run using parameters that caused it to exhibit 
m—~~cmaraing force corresponding to the curve at sea state 
four, therefore this curve was used to remodel the RTSSD. 
The new simplified equation of motion for surge acceleration 


became 


me fn) - (FSD/m) ev. +r 


forw 


Mmere FSD 1s full scale drag which 1s defined by the curve 
moierre. 5. TEORW is the sumation of the four thrusts assumed 
in this model to be applied at the same point along the 
eemcer line of the craft. 
2. Implementation 
FSD was implemented into the RTSSD model by dividing 
the curve into four separate regions and defining FSD within 


each region. In region I the retarding force was approxi- 


mated by a normalized velocity cubed term. In regions II 
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and III straight line approximations were used. In region 
IV a polynomial curve fit was used. The equations defined for 


each region are listed below. 


Region I - ( 0-18 knots) 
ie 3 
ee CDx4 ( u/ Vmax ) 
Region II - (18-22 knots) 


FSD 


CONSTANT, ~ “pxs u 


Region III - (22-36 knots) 
FSD = Cnxe u - CONSTANT 


Region IV - (above 36 knots) 


- 2 ile) 


Z 


Seeelaentitication of Surge Drag Coefficients 


iheregion 1, °C is equal to the first peak value 


Dx 


miecanag. In regions II and III, C and Cox? are ee dude 


DXS 
the slopes of the drag curve and the constants were solved 
Pmepy appiying the boundary conditions at the intersections 
Beene adjacent regions. The values of Chx3 and Cyox4 were 


computed by taking FSD at two different known values of u 


and solving for Coys and Cox Simultaneously. 


Pee LENEARITZED MODEL ANALYSIS 
Simce the surge model of the RTS5D is extremely non- 
linear over the complete speed range, analysis is very com- 


plex. Therefore the model was linearized and it's performance 


38 








analyzed at various points on the FSD curve. Linearization 
of the model was accomplished by computing the linear 
coefficients through partial differention with respect to 
surge in each of the four region equations and evaluating 
these terms at various operating point surge velocities (u.). 
The Linearized model is shown in Fig. 6 and the linear coef- 


ficients for the surge drag force FSD are shown below. 


Region 
: PSD ( ) = 3 Cyy, ( uo/ a 
ta FSD (4 ) a igs 
rt FSD, ‘ 


O 


/ 


O 


: 2.5 
2 Cpx3 Yo - 145 Spy / Uy 


A map of the linearized surge directional roots are 
Bown in Fig. 7. The characteristic roots for each operating 
region are shown coded (see legend). Note that only Region I 
and Region IV roots are u, dependent and therefore move as US 
changes. The direction of change along the real axis ais 
noted by the arrows shown. Region II and III roots are fixed 


in locations and independent of u.- 
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Figure 6 


Linearized Surge Model 


Region [IV of Fig. 7 is of particular interest since the 
ability to operate in this speed range is a primary consi- 
deration in the development of the SES. It can be seen 

from the map of the roots that within this region, at speeds 
from 36 to almost 60 knots, poles exist in the RHP. This 

is a result of the negative drag slope where the retarding 
force actually decreases with an increase in speed, the 
magnitude of variation being dependent on sea state. Region 
Mmimauso exhibits this negative drag characteristic but since 
7€ is Of such a narrow speed range, operation in this re- 
gion Was not considered in this design. Further modeling 
changes for the other degrees of freedom were based on 

Region IV response. In Regions I, III and in Region IV above 
60 knots, the system will be stable at the speed US where FSD 


1s equal to commanded thrust. 
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Figure 7 
Linearized Surge Root Map 
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[i it were necessary to operate at a speed within the 
negative drag slope characteristic, it would be necessary to 
actively modulate the thrust command in order to maintain 
a nearly constant speed, but this would require a great deal 
of operator attention since the negative drag slope in Region 
IV is one of slowly divergert instability. For this reason 
a closed loop speed controller has been included in the data 


base model and therefore was also added to the RTSSD model. 


foo IMPLI FIED PROPULSION DYNAMICS 

Before adding a speed autopilot to the model one addi- 
tional refinement was included in the design. Thrust command 
changes in the RTSSD as reported in Ref. I were instantan- 
eously entered into the equations of motion as thrust 
Changes. This was not the case for the data base simulation 
program. Specific thrust levels are generated by specific 
water flow rates through water jet nozzles. These flow 
rates are nearly proportional to the pump speed that pro- 
duces them. Since the pumps are driven directly by gas 
turbine engines, a change in thrust occurs with a change in 
turbine speed. <A thrust command translates to a power tur- 
bine speed command which as an output has power turbine 
speed achieved. This output is translated into thrust 
achieved. It was specified in Ref. 7 that one of the goals 
in the base line design was to have power turbine speed 
Fomtro, response charasteristics which are relatively inde- 


Pendent of the power setting. Accordingly, in the design 
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of the auto pilot of Ref. Pee ie simp lit iod souoomes aon 
D lock shown in Fig. 8 was used. The magnitude of the pro- 
pulsion system time constant, t was not specified. For 


the purpose of this analysis +t was chosen to be 1 second. 





Figure 8 


Simplified Propulsion System 


Wee ier SiS OF LINEARTZED MODEL WITH SPEED CONTROLLER 

The complete linearized model with speed controller that 
Bemecluaded in the model for the purposes of this study is 
Shown in Fig. 9. The saturation effects of fuel control are 
ignored and the linear representation is shown as Ke which 
is assumed to be unity. The controller includes propor- 
mlotialeplus integral control. The proportional control was 
included to keep the ship speed error low while the integral 


“OienOlewas included in order to zero out long term control 


caer Or . 
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The characteristic equation for the system in Fig. 9 is 
3 : 2 ; =) 
S- + (1+ FSD), )/m)S + 1/m (FSD), ) + K) S+ K/ mf 
0 0 


Using the Routh Hurwitz Criterion for stability 


53 1 1/m (FSD(y ae ) 
a (1+ FSD (4, y/ m ) K/ mT 

st 4/m (1 + FSD (4 y/m) (FSD(, er) Kem te 0 

30 Ket Fe 0 


For Stability K must be greater than O and greater than 


42 } 
PSP (u.) T + m7 FSD /,, ) 


O 


/ 


ime Values of controller parameters K and t can be initialized 
by the operator for either loosely controlled speed or tightly 
controlled speed for such situations as station keeping, 
PereEPs, Or test operations. For comparison to the data base, 
these parameters were selected to achieve closest agreement 
between the values of surge achieved by the RTS5D and DBSIMSD 
models after the completion of a 360° turn with 50 knot initial 


Surge veloctiy and a thrust effector angle of 15°. This test 
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condition for the nonlinear models was chosen for two reasons. 
First, given fuel conservation as a necessity, 50 knots is a 
merely heh speed initial condition because it is close to 
minimum drag. Second, this initial condition provided an 
opportunity to test the operation of the speed controller in 
the negative drag region of FSD since the decrease in speed 
going into the turn with u,=50 knots caused increased drag. 
The values of controller parameters that caused closest agree- 
ment between the two nonlinear models, when tested in the 
linearized model showed instability in the linear sense. 
This is desirable in a turning mode for the nonlinear model 
Decause it allows speed to decrease and the turn to be com- 
meeetcammoere quickly. Stiff speed control could be achieved 
by increasing the controller gains, but would normally be 
used only in straight ahead runs. 

Coupling of the equations of motion developed in Chapter 
II necessitated a repetition of the parameter identification 
Lue ande-@ . 


Dy DP DAP 


These values were adjusted for the best curve fit compared 


Procedure resulting in new values for C 


to the data base using the linearized model for the RTSSD 
moe oo0 turn using 15° thrust effector angle with initial 
surge velocity 50 knots. An identical series of tests to 
those in Chapter III were conducted comparing the new model 
of the RTSSD to the DBSIMSD. Complete validation results 


are tabulated in Chapter VI. 
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Pee crECTIVE THRUST EFFECTOR ANGLE 

As a result of the negative drag characteristics imposed 
on the model, RTSSD yaw rate (r) was high in the low speed 
test and low in the high speed test. This opposite error 
at both ends of the test speed spectrum was reduced by 
introducing a correction factor to the thrust effector angle 


proportional to yaw rate 


2 off =e ee SLIP = rx 


and adjusting 2£w to increase r. This had the desired effect 
of reducing error at both ends of the test spectrum since 
the reduction of the thrust effector angle was greater at 


the high end. 


pee ekOACH CONDITION FLAG 

As stated in Chapter I, broach response characteristics 
@eeemoe included in this model, therefore a broach condition 
warning flag was included in the design in order to alert 
the operator who could then take corrective action to avoid 
broaching. 

Broaching is caused by a water intake being exposed to 
air and is a function of surge speed, roll angle, pitch 
mete, sarint angle and sea state. For the data base 3K=SES5 
with broach boundaries defined at 10% air by volume, the 
Deoatchaumeas boundaries in calm water in terms of roll angle, 


peer angle and drift angle are shown in Fig. ll for 40 knots 
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moet. 12 for 60 knots. Reference 6 listed maximum thrust 
effector angles as a function of surge speed and air flow 
mieerder to avoid broaching. This condition is shown plotted 


in Fig. 10 for fixed plenum air flow rate. 


EFFECTOR 
ANGLE 
(DEG ) 


20 


MS 
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BROACH 
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Figure 10 
Maximum Thrust Effector Angle As A Function 


Of Surge Velocity In Order To Avoid Broaching 


As stated in Chapter III the model developed in Ref. 1 did 
not flag the broaching condition. To implement the broach 
warning flag and display the other new parameters for the 
modified RTS5D the pilot graphic display described in Ref. 1 
was changed to incorporate the necessary data. Figure 15A 
depicts typical information displayed to the pilot under 


normal operating conditions. The new information to the 
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ieraeer siown in the lower right of Fig. 13A 1s sea state 
foe ol), Computer assist (CASST, on=1 off=0), controller 
gains (K,KK) and speed commanded (SPCOM). [If the operator 
Eeaeeccacmtene thrust effector limits (ZMAX) to prevent 
broaching (see Fig. 10), a flashing warning is displayed 
mimabace Of the bottom Line of information on the screen. 


iis) is shown in Fig. 13B. 
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VARIABLE DYNAMICS SEA TRACK 


U (HEAVY) 
PpOl LIGHT) <= 


V (HEAVY) 
VDOT (LIGHT) [N= 


R (HEAVY) 
RDOT (LIGHT) 





P (HEAVY) |. 
PpOT (LIGHT) F® | 

| 
Q@ (HEAVY) |. _ 1 
QDOT(LIGHT) : 

\ 

CONTROL INPUTS TIME NAVIGATIONAL DATA 
THRUST RUDDER DELT ELAPSED PITCH HEADING/RATE DRIFT ROLL SPEED 
300000 -.20 .084 eee , 52 325/2.1 -,28 8 “S@Pi 
[7 78 T9 T10 SEAST CASST K KK  SPCOM 
76000 75000 75000 750000 4 1 2900 3 50 


Figure 13A 
Pilot Graphic Display (Normal) 
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VARIABLE DYNAMICS SEA TRACK 


VDOT (LIGHT) &e 
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R (HEAVY) WA ee ee ae es 
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ie (HEAVY) 
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i‘ 

CONTROL INPUTS TIME NAVIGATIONAL DATA 
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—— 2 Oo fC Hn 

/ ae 

Figure 13B 


Pilot Graphic Display (Broach Condition Exists) 


52 





V. IMPLEMENTATION 


A. INTRODUCTION 

Given the existence of a Real. Time Simulation of a SES, 

a decision was made to refine that model using more accurate 
equations of motion, improved computer iteration time and 
operator hardware modifications. A complete description of 
mace RISO) is contained in Ref. l, therefore only changes 

are discussed here, and included are the introduction of 
negative drag in a sea state, propulsion dynamics, speed 
Someroaleand a broach condition warning in software. Hard- 
Ware modification included changing the thruster control 


box to make it more accessable to an operator. 


B. REQUIREMENTS 

The following cra@teria were required of the remodeled 
aco Le 

1) The refined simplified equations of motion would be 
solved and the results output on a real time basis. 

2} The solution would be subject to real time control 
efforts generated by an operator observing the output. 

3) Computer iteratiom time was required to be less than 
100 ms. 
A more complex discussion of condition 3) is presented in 


Appendix A. 
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C. HARDWARE DESCRIPTION 

The hardware arrangement. in block diagram form is shown 
in Fig. 14. The only change to the RTSSD was the redesign 
emo relocation of the thruster console to the pilots seat 
mere sit Can be more conveniently utilized by the operator. 
The thruster console uses "linear movement." potentiometers 
in place of previously used "dial" potentiometers in order 
Bomnonre realistically simulate thrust controls for the SES. 


fjewcnruster console is shown im Fig. 15. 


D. SOFTWARE DESCRIPTION 

The remodeled software package was a FORTRAN IV digital 
program which consisted of a main graphics program and three 
major subroutines. It was necessary tio make a third sub- 
routine out of the equations of motion computation that. had 
previously resided in the main program prior to making any 
Changes to these equations. It was discovered when the 
first program modifications were attempted that the compiler 
for the XDS-9500 computer was operating at its limit with 
the program as it then existed, therefore, before changes to 
the equations could be implemented, the program structure 
had to be modified. The failure mode was eliminated from 
the new model because of initial assumptions that the thrusts 
would be all applied at the same point. Thrust failures 
Should not be considered until the model includes broaching 


dynamics which is a major source of thrust loss. The flow 
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chart for the remodeled RTSSD is shown in Fig. 164. The com- 


plete multiplexing algorithm is shown in a flow chart in 


Fig. 17. 
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Figure 14 
RTS5D Mod Block Diagram 
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Figure 15 


Thruster Console 
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INITIALIZE 
INPUT 
“ CHANGES 


> DARD SecmceKs 


1. CALL SUBROUTINE AD 
2. CALL SUBROUTINE ONEMOT 
3. CALL SUBROUTINE BOAT AND DISPLAY 


GRAPHICS DISPLAY MULTIPLEX ‘ALGORITHM (FIG 17) 
1. READ CLOCK SSE? STEP SIZE 
e. RESEEeCLOCK 
C8 lao > LE af 
=f} “ye 


NO 


Figure 16 
RTS5D Mod Program Flow Chart 
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INCREMENT MULTIPLEX GRAPHICS DISPLAY INDEX I 


DISPLAY INFO (I) ON 
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ENCODE GRAPHICS 
INFO (I) 
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Figure 17 
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RTS5D MOD DISPLAY MULTIPLEX ALGORITHM 
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VI. kYSSD MODS I AND II RESPONSE CHARACTERISTICS 


RTSSD MOD I includes negative drag effects in a sea state, 
simplified propulsion dynamics, speed control and a broach 
condition warning flag. RTSSD MOD II is identical to MOD I 
Saeepe Gilat it includes the correction factor for thrust 
effector angle developed in Chapter IV, Section E. 

The validation of the modified RTS5SD was accomplished 
in the same manner as the RTSSD. (See Chapter III). The 
results are displayed in the following tables. <A complete 
560° turn comparison of the DBSIMS5D, the RTSSD and the modi- 
meade ktoosD at 50 knots, 15° thruster angle in shown in Fig. 
ig. Additionally, a comparison of transient response was 
undertaken and the time of first peak for the variable v, r, 
and 6 for each model was graphed in Figs. 19-21 for speeds 
Goes and 60 knots. Transient response characteristics 
are an important consideration in the development of a real 
time model and are discussed in Appendix A. The transient 
response of the RTSSD models was faster than that of the 
data base model with the exception of roll angle (¢#). This 
mmeacted on the validity of a real time solution of the equa- 


Biohicmor motion and the display of the results. 
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TABIE™ TV 


FesoobD Mod I and Mod IE Performance Test? at 40 Knots 


Bac ector 
angle RTSSD Mod I RTSSD Mod II 
> in degrees $ 
ist p.K Qss Ist pk Qss 
Se {2 ~45 . 40 oO 
- . 96 92 .79 .76 
8 imac ces 6 
esses p< Qss isco Qss 
5 ee 2 lela i 20 iS 
0 20 ies eZ elke 
u ive G/ Sec u 
ste pik Qss FS tpi Qss 
ae ee DO 62704 62.58 Oy 05 
EO) ° n/a leo! n/a Oy 
V ine, SEC V 
Sic pik Qss Aston (1K Qss 
5 2, (8) 7s 2d Zo DT 
0 eo Sierore So Sa.2 
ii in degrees/sec fr 
lst pk Qss Sete p< Qss 
5 .69 59 sore ao 
1 Oe eS de 2 1.04 1.09 
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TABLE V 


PieoooUeMed [T and Med II Performance Test at 50 Knots 


effector 
angle RTSSD Mod I RTSSD Model 
> in degrees $ 
lst pk Qss Sp Qss 
Sig Bei) 39 aS mSo 
O° . 90 oe . 85 nS 
eS? 5) £5.30 1.24 1.19 
6 in degrees 6 
lst pk Qss ES pi Qss 
oe eee 1a 9 eS a2 0 
FO'° eaZ4 egAleS 1.25 tec 0 
15° 1.24 ene es L219 
u in ft/sec u 
Est) pk Qss lst pk Qss 
Se Ty a Ses mya 82.98 
O° n/a 78.68 n/a 79.55 
eS n/a ooo n/a Leas 
Vv if tty see V 
lst pk Qss sy p< Qss 
Se Do, 250 2.39 De 9 
10° 3.48 33 See 5 Sl 
5S” aS 29 3.98 4.00 
ve in degrees/sec ye 
ee) pik Qss lst pk Qss 
SS no 39 49 41 
Oe 96 o5 . 88 oe 
eS be S ey Heese 4 





TABLE VI 


ploob Mod IT and Mod II Performance Test at®60 Knots 


efttector 
angle Rissaiode. RTSSD Mod II 
b in degrees > 
lst pk Qss lst pk Qss 
5 LoS ao 46 v4 
8 in degrees 6 
lst pk Qss lst pk Qss 
5 eee 5 ieee ea ec 
u u 
lst pk Qss lst pk Qss 
Ss m/ a 10 Umea my a 99.59 
V V 
is Gap K Qss its Gena Qss 
5 2.64 aoe 2.49 LoS 
i i 
lst pk Qss lst pk Qss 
5 oie oe 46 O06 
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Wie 50 knots 
15 (effector angle) 
10 second intervals 


RTS5D 
DBSIM5D 
RTS5D MOD I 
RTS5D MOD II 


27> 


x * BO 


Figure 18 


360 Turn Comparison of Four Models 
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BS 40 knots Q@= RTSSD 
fffector angle = 5 e Eee a 
Figure 19 


Response Time Comparison at 40 Knots 
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as 50 knots ©= RTSSD 
Effector angle = 5 — AE i 
Figure 20 


Response Time Comparison at 50 Knots 
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Figure 21 
Response Time Comparison at 60 Knots 
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Vi eee ON CIEISEONS 


The RTSSD has been shown in Section VI to be a viable 
input/output model of a 3K-SES by closely matching its output 
Snaracteristics for a given input to those of the data base 
model. 

iiemrerinement of the right side of the simplified equa- 
tions of motion, especially the inclusion of the negative 
drag characteristic in a sea state has made the model more 
accurate over a wide speed range and gives the operator, upon 
program initialization, the ability to choose different drag 
Characteristics caused by sea state. 

The model was tested in a negative drag region with 
automatic speed control amd simplified propulsion dynamics 
and it was demonstrated that by properly adjusting controller 
gains that the output characteristics of the DBSIMOD were 
Closely matched. 

Previously conducted thrust failure analysis with man 
generated corrective action has been shown to be invalid 
Until broach dynamics are included into the model. Previous 
tests allowed maximum thrust effector angle correction regard- 
less of this consideration. The modified version of the RTSSD 
does not include complete broach dynamics but does warn the 
operator if he should exceed a maximum thrust Ee eetor anne 


for a given forward speed. 
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The 3K-SES may be operated in the negative drag slope 
Speed region approximately (36-58) knots without the aid of 
a speed controller but would require constant manual modula- 
tion of the thrusters. Since minimum manning is one of the 
criteria in 3K-SES desigt, a speed controller is necessary 
to free operator attention to other duties when operating 
in this speed range. 

Another consideration of the implementation of the non- 
Linear drag curves is that for a given sea state, two values 
of minimum drag can be determined, one in Region IV ard one 
between Region II and Region III. Knowledge that drag can 
be lower in a particular operating region for fuel conserva- 
tion is of prime importance. 

It 1s recognized that a time delay exists in the solution 
of the equations of motion and a time delay exists in the 
eaapntcal display of the solutions to the operator. ReaL 
time interaction between the man in the loop and the computer 
is achieved by setting the integration step size equal to 
the computer loop iteration time. Real. time considerations 


are further discussed in Appendix A. 
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VIII. RECOMMENDATIONS 


The following recommendations are made concerning future 
improvements in the RTS5D model. 

1. Further refinement of the right hand side of the 
equations of motion is necessary in order to more closely 
match the characteristics of the data Base model. 

2. Implement software broach dynamics in the model with 
thrust applied differentially at different lever arms Csee 
Fig. 3) so that thrust failures and effectiveness of correc- 
tive action may be evaluated. 

3. Rewrite the program such that the AGT-10's would 
provide asyncronous parallel graphics Cexisting equipment) 
or implement the model on a more modern system which uses 
parallel processing and improved graphics. 

4, Introduce the sixth degree of freedom Cheave) as a 
parallel operation. Since heave time constant is much faster 
than the rest of the system dynamics, real time heave analysis 
could be conducted simultaneously if cross coupling is neg- 
wected. 

9. The system is designed for man in the loop interface, 
so further operator testing is recommended as a valuable 


SoOuree Of information for future refinement of the model. 
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APPENDIX A 


REAL Pee oy Slo 


i eos Nelson, in Ref. 1, accomplished real time solution 
6f SES equations of motion with the use of a 1000 Hz clock 
to measure computer iteration time and by using that time as 
an integration step size for rectangular integration to solve 
the simplified dynamic equations. He concluded, by comparing 
th. solutions of these equations to that obtained by using a 
Runga Kutta variable integration step size, that the computer 
meerac1on loop time must stay below 100 ms in order to main- 
tain at least 3% accuracy of solution. He further developed 
a multiplexing algorithm in order to output these solutions 
tO an operator in the form of graphics on a real time basis. 
This model can be invalidated as a real time model in two 
ways. First, in order to have real time, the integration 
step size used must be equal to the computer iteration time 
and in this model, because of the multiplexing algorithm, 
each loop time is different. The integration step size, 
dt determined by measuring loop N is used to solve the equa- 
tions in loop Ntl, which has a different loop time, there- 
fore the solutions are not exactly real time. Secondly, 
the multiplexing scheme makes the model into a sampled data 
system with sample time (Ts) approximately Z seconds. Even 
inane first condition did not apply and real time solutions 


to the equations could be assumed, then the output to the 
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operator would be sampled real time data. Given the fastest 
mimeeconstant for the RTS5D simplified equations (approxi- 
meeely 1:3 seconds for pitch rotational motion), it can be 
Seemethat the Nyquist Criteria is clearly violated. Yet 
even with these obvious limitations, it has been demonstrated 
that the RTSSD can closely reproduce the output characteris- 
tics of the DBSIMS5D and therefore is a useful model. Surge 
speed as well as sway, yaw and roll (Figs. 21-23) have much 
longer effective time constants. It is recognized that 
there is a delay between the real time of any event and the 
observation of that event, therefore "Real Time" as used in 
this report is a matter of degree, or more specifically 

how much delay time between solution and display can be 
tolerated. Given the scale of Praanaas presented to the 
operator in this model and the speed with which the SES 
responds, the RTSSD is an acceptable model, however it 

would be desirable to reduce the sample time ae and make the 
integration step size equal to the computer loop iteration 
time for each iteration. 

An algorithm to set the integration step size equal to 
the computer iteration time of the current loop was imple- 
mented as a software design. cdrange to the RTSSD model. This 
was accomplished by defining the integration step size 
as an array DELTCN) with N equal to the total number of itera- 
tions nNe.cessary-to display all of the graphical information 


to the operator. Each loop iteration time was.m.asured and 
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stored in this array during the first pass through the 
multiplexing algorithm. DELICN) was initialized at zero 

so that time would not increment and integration would not 
be performed during the first pass through the multiplexing 
algorithm. On subsequent passes through the multiplexing 
algorithm the integration step sizé was equal to the loop 
Hime in which it was used. 

This timing algorithm wren tested in the RTS5SD model 
was found to accurately represent real time but the devia- 
tion from previously used integration step size values 
Seaiaced instability 1n bothwthe pitch and roll dynamics. 

This needs further investigation. 

Sturgeon, in Ref. 8, demonstrated the valw. of using 
poralilelwerray processors to reduce computational time in 
real time simulations. Such a model could be designed for 
the RTSSD using the existing equipment. The AGT-10's which 
are. used for displaying the output can interface with the 
XDS-9300 main computer in two ways. The multiplexing scheme 
of the RTSSD uses a "hand shaking" method with a GRAPHO or 
TEXTO subroutine call in each loop. It 1s precisely the 
inclusion of graphics in each loop that causes problems for 
real time presentation. The solution of the equations requires 
only approximately 35 ms. Each GRAPHO or TEXTO call, with 
its associated "hand shaking" time, makes up the rest of the 
loop time, approximately 40 ms more. It can be seen that By 
introducing two such calls in one loop that the 100 ms time 


restriction would be exceeded. 
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An improvement to this method would be to use the second 
interface between the AGT-19's and the XDS-9300 which is 
thet the AGT-10's have a direct connection to some memory 
locations in XDS-9300. By solving the dynamics for the SES 
and storing the data in these memory addresses in the form 
of arrays, the AGT-10's could ke. made to operate asyncronously 
in parallel as array processors by using a process known as 
"cycle stealing” to fetch the data from the XDS-9300 and to 
display this on the screen for the operator. This process 
would have the advantage of reducing and making more uniform 
the computer iteration time used as the step size for the 
rectangular integration in the solution of the SES equations 
of motion. The sample time would also be reduced because all 
of the time wasted in the “hand shaking” is eliminated. It's 
drawback is that the output is not truely real time because 
of the asyncronous operation of the AGT-10's. However, the 
reduced sampled data time and computer iteration time would 
make it a closer approximation of real time than presently 


exists. 
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pe Ls 
ARCLL 
ASTOP 
AXST 
Eel L 
esta 


BOAT 


APE aN OMeN 5 


RTS5D MODIFIED PROGRAM NOMENCLATURE 


A/D trunk 500 voltage 

A/D trunk 501 voltage 

Added mass coefficient in yaw moment equation 
A/D trunk 502 voltage 

Added mass coefficient in roll moment. equation 
Added mass coefficient in pitch moment equation 


4/D trunk 503 voltage 


array of A/D lines 

speed line index counter 

mass 

pitch angle in degrees 

4 iterations between print execute commands 
Real IARR (1) 

turn rate in degrees 

roll angle 

performance index J, Limit 

initial condition X-coordinate on sea track 
igi ral condition Y-coordinate on sea track 
demre. angele 


subroutine to generate perspective of SES 
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e2aSst 
Spr 
(4) 
GY 
SDZP 
DAL 

Be LT 
DELTA 
DFCCM 
DFCCM 1 


DFCCM 2 


eriMe 
BOIMPLT 
DUCCM 
BALLOT 
meee LOT 
Deen er es 1. 
EFLG 


Bie We 


automatic speed control mode indicator 

dow real forces lumped coefficient 

drag forces lumped coefficients 

sway forces lumped coefficient 

Sidewall rolling moment lumped ceofficient 
iia! COeamaloge call 

iteration loop time 

commanded iteration loop time 

total change in fuel command from speed comtroller 
change in fuel command due to proportional speed 
controller 

change in fuel command due to proportional speed 
controller 

Vieaeaeme 

scaled value for v 

difference between ucom and u 

same as DTIMPLT 

scaled value for v 

restart value for v 

emmoreeriag for TEATOC call 

r time 

restart value for r 

scaled r time 

scaled r 

stern buayancy force 

bow seal pitch force 


bow bduoyance force 
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POIME 
FTIMPLT 
FXPLOT 
FYPLOT 
FYREPET 
é 
GRAPHO 
GTIME 
GTIMPLT 
GXPLOT 
PLOT 
GYREPET 
Gy 
HEAD 
HOLD 
HTINE 
HTIMPLT 
HYREPET 
HYSESX 
HYSESY 
TARR(1 ) 
FeAaeeri( 2 ) 
TARR(3) 


fuel command 

port sidewall buoyancy force 
starboard sidewall buoyancy force 

v time 

scaled v time 

same as FIIMPLT 

scaled v 

mecwarteyaiue for v 

gravitational acceleration 
subroutine to project image (non alpha-numeric) 
R time 

scaled R time 

same asGTIMPLT 

scaled R 

restart R value 

dynamics maximum and minimum array 
craft heading 

subroutine to freeze display/program 
u time 

scaled u time 

analog restart y coordinate on sea track 
Anatomy position of craft 

analog x position of craft 

digitized mode selector 

digitized effector angle 


GQrewtt Zed. aa 


we 





TARR (4) digitized 78 


TARR(5) Geer zed C9 

fen (6) digitized 1T10 

TARR(7) digitized surge velocity command 
Lee, multiplex graphics index 
ICOA speed line index 

moRAFT vehicle designator 

TER exmror ilag 

PPR memes history Lama 

fee RD time history limit 

EIPRE tame history limit 

TIPRE mene Hescory Limit 

LIPRG Time wiwvstoryelimat 

ie RY tame history limit 

Hee RC time history limit 

ire R P mamee 1 Story limit 

ier Q time history limit 

PEP RR time history limit 

PRS Timer nl story lim t 

Teme Binienc COUR Ter 

Ti craft perspective counter 
IJA craft perspective counter 
IJAD time history counter 

LJAE Pines nr Stony counver 

IJAFP time history counter 

ITJAG time Mieeore counter 

I JAH time history counter 
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IJAO time history counter 
IJAP time history counter 
IJAQ time history counter 
IJAR time history counter 
IJAS time history counter 
IJD time history counter 
IJE time history counter 
IJF time history counter 
IJG time history counter 
IJH time history counter 
IJO time history counter 
IJP time history counter 
IJQ time history counter 
IJR time history counter 
IJS time history counter 
IL graphics console display array 
ILAL graphics console display array 
ILA2 graphics console display array 
ILA3 graphics console display array 
ILA4 graphics console display array 
ILA5 graphics console display array 
ILA6 graphics console display array 
ILA7 graphics console display array 
ILA8 graphics console display array 
ILA9 graphics console display array 
ILA1O graphics console display array 
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eA 1 
iA 1 2 
fines 
ILA14 
TLAIS 
ILAI6 
ILAL7 
Ives: 
ILAIS 
2 0 
ILA21 
ILAZ2 
ILA23 
ILA24 
TLA25 
TLA25 
Th A27 


IL A28 


Aco 


IL AP 
IPACK 
mee OT 
Eee 


IPRD 


graphics 
graphics 
graphics 
graphics 
graphics 
graphics 
graphics 
graphics 
graphics 
graphics 
graphics 
graphics 
graphics 
graphics 
graphics 
graphics 
graphics 
graphics 
graphics 


graphics 


subroutine to load display arrays 


console display 


console display 
console display 
console display 
console display 
console display 
console display 
console display 
console display 
console display 
console display 
console display 
console display 
console dispLay 
console display 
console display 
console display 
eemsole display 
console display 


console display 


array 


array 


array 


array 


array 


array 


array 


array 


array 


array 


array 


array 


array 


array 


array 


array 


array 


array 


array 


array 


graphies console display array 


Zilomimts hom aecotnner LIM t 


time history counter limit 
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PERE 
Pe 
IPRG 
PPR 
ee RO 
ie RP 
EP RO 
UIE ase 
eS 
ISWA 
Pe 
IVEIWA 
IVIEWAA 
IVIEWB 
IVIEWC 
IVIEWD 
IVIEWE 
ny LEEW 
IVIEWG 
IVIEWH 
IVIEWO 
IVIEWP 
IVIEWQ 
IVIEWR 
IVIEWS 


IVIEWZ 


time 
time 
time 
time 
time 
time 
time 
time 


time 


display array 
Subroutine to display alpha-numeric data 


craft perspective/north road array 


ha Stony 
AUSEOLY. 
Rustony 
history 
history 
history 
Piston, 
Ras cain 


history 


counter 


counter 


counter 


Counter: 


counter 


counter 


counter 


Counter 


counter 


sea track array 


Sea track axis 


array 


u time 
v time 
r time 
v time 
r time 
u time 
D time 
q time 
p time 
q time 


analog 


Ais t@isy. 
histormy 
history 
history 
history 
history 
history 
history 
history 


history 


array 
array 
array 
array 
array 
array 
array 
array 
array 


array 


ila jaleiys 
earn at 
ee 
Ilse 
eerie te 
Pama 
Ihotga(a bie 
Logik 


eterna 


sea track array 


time history axis array 


Se 





NAD 
ONEMOT 


Mane oe momen: Clo Inert ula 


time 
Tene 


Wale 


time 
time 
time 
time 
time 


time 


ALSvory 
fis voOru 
ies Om 
histroy 
history 
history 
ha srocy 
history 
history 
Rene aor, 


nestor 


Coun ves 
counter 
counter 
counter 
counter 
Councer 
counter 
counter 
counter 
Counter 


counter 


Peepoeruronal speed controller gain 
intregal speed controller gain 
lever arm for buoyancy force 

draft of sidewall 

Side thrust lever arm 

lever arm for sway drag force 

Cloece Inverrsupe Count 

AD trunk line array 

Sueseutine wo Calculate SSeS dynamics 
lever arm for thrust side component 
5 time 

scaled D time 

same as OTIMPLT 

scaled D 


restart p value 
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'J 


ko 


QDcT 
QDCTH 
QMAX 
QTINE 
S2iNPLT 
QAPLCT 
QYPLOT 
eer P eT 
R 

RDOT 
RDCTM 
READCLOCK 


roll rate 

plenum pressure 
BoOliwacee Leration 
Maer oll accsleration 
roll angle 

max roll rate 
craft heading 

q time 

scaled q time 
same as PTIIMPLT 
scaled q 

mes Cai q value 
Peper rate 
Wiesemuacce |eration 
max q 

max q 

p time 

scaled p time 
same as QTINPLT 
scaled p value 
restart p value 
ena luna rave r 
r 

max r 


subroutine to sample real time clock 
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“Ssyoak 
B10 
RMAX 
aol i 
oliver LT 
RX 

mA LCT 
cog 

aoe LCT 
monks 
mz 


SPDLIN 


oe 


STARTCLCCK 


SLIME 
SlLMPLT 
SAPLCT 
sa PLCT 
Sere T 
Sian PET 


subroutine to reset analog computer 
density of water 

maximum yaw rate of craft 

q rate 

scaled q time 

ly 

same as RTIMPLT 


ly 


— 


scaled gq value 
restart q value 
=a 
sea state 

scaled we 

scaled ty 

scale factor in time history arrays 
scale factor in time history arrays 
Beta influence coefficient 
sceedline lower screen limit 
Gpanets Velocity in knots 
subroutine to start clock 

analog time 

scaled analog time 

EXYSESA 

AVS one 

restart value for HYSES&’ 


Gestalt Value «.Or Histol 
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210 
TLOMAX 
at 
T7TAS 
we 
T8MAX 
T'9 
TOMAX 
age 


eDCT 


TIE 

Pp LT 
INT 
ee iLe 
eT LEO 
ere 1 
TM 


ECetl 


total thrust 

thruster number 4 

max TiO value 

thruster number 1 

max T7 value 

thruster number 2 

max value of T8 

Gaeuste> number 3 

max value of T9 

subroutine to display alpha-numeric data 
T 

pitch angle 

u time 

time 

sealed TP 

real clock interrupt frequency 
program name display array 
program name display array 
program name display array 
manual thrust command 

max T 

sum of side thrust components 
sum of yaw moment 

forward velocity (surge) 


forward velocity (surge) command 
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eRe PET 
foot 
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u 

max u 

max u 

scale factor in speedline generation 
side velocity (sway) 

thruster console pot array 

“ 

max v 

max v 

total velocity 

heave velocity 

heave acceleration 

width of dow seal 

subroutine to assign clock interrupts to 


variable 


lever arm for sway drag forces 


boat verspective scale factor 
boat perspective scale factor 
display array 

reset value for SESX 


Seale factor for analog plot 
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Appendix D 
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